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^ ! Abstract 

We study the evolution of shear-free spherically symmetric charged fluids in gen- 
eral relativity. We find a new parametric class of solutions to the Einstein-Maxwell 
system of field equations. Our charged results are a generalisation of earlier treat- 
ments for neutral relativistic fluids. We regain the first integrals found previously 
for uncharged matter as a special case. In addition an explicit first integral is found 
which is necessarily charged. 
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1 Introduction 



x. 

Solutions of the Einstein-Maxwell system of equations are important in rela- 
tivistic astrophysics as they may be used to describe charged compact objects 
with strong gravitational fields such as dense neutron stars. Several recent 
treatments, including the works of Ivanov p] and Sharma et al [2\, demonstrate 
that the presence of the electromagnetic field affects the values of redshifts, 
luminosities and maximum mass of a compact relativistic star. The electro- 
magnetic field cannot be ignored when considering the gravitational evolution 
of stars composed of quark matter as pointed out by Mak and Harko |3J and 
Komathiraj and Maharaj [3]. Therefore exact models describing the forma- 
tion and evolution of charged stellar objects, within the context of full general 
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relativity, are necessary. Electromagnetic fields play a role in gravitational col- 
lapse, the formation of naked singularities, and the collapse of charged shells of 
matter onto existing black holes (as indicated by Lasky and Lun 05]). Signif- 
icant electric fields are also present in phases of intense dynamical activity, in 
collapsing configurations, with time scales of the order of the hydrostatic time 
scale for which the usual stable equilibrium configurations assumptions are not 
reliable (as shown in the treatments of Di Prisco et al j7] and Herrera et al [8]). 
It is interesting to note that Maxwell's equations play a role in several other 
scenarios, including the evolution of cosmological models in higher dimensions. 
De Felice and Ringeval [5] considered braneworld models, exhibiting Poincare 
symmetry in extra-dimensions, which admit wormhole configurations. 

Spherical symmetry and a shear-free matter distribution are simplifying as- 
sumptions usually made when seeking exact solutions to the Einstein field 
equations with neutral matter. The field equations may then be reduced to 
a single partial differential equation. What is interesting about this equation 
is that it can be treated as an ordinary differential equation. A general class 
of solutions was first found by Kustaanheimo and Qvist [10]. Comprehen- 
sive treatments of the uncharged case are provided by Srivastava [11J and 
Sussman [12]. The generalisation to include the electromagnetic field is eas- 
ily performed and is described by the Einstein-Maxwell system. The field 
equations are again reducible to a single partial differential equation, now 
containing a term corresponding to charge. A review of known charged solu- 
tions, admitting a Friedmann limit, is contained in the treatment of Krasinski 
|13j . A detailed investigation of the mathematical and physical features of the 
Einstein-Maxwell system has been performed by Srivastava [H] and Sussman 
[T5|[T6] respectively. 




The objective of this paper is to investigate the integrability properties of 
the governing partial differential equation that contains a term correspond- 
ing to charge, for shear-free fluids. This investigation is performed using an 
elementary approach suggested by Srivastava [IT] . In Section 2, we reduce 
the Einstein-Maxwell field equations, generalising the transformation due to 
Faulkes [17J, to a single nonlinear second order partial differential equation 
that governs the behaviour of charged fluids. As in the uncharged case, this 
equation can be treated as an ordinary differential equation. In Section 3, we 
derive a first integral of the governing equation by generalising the technique 
of Srivastava [UJ first used for uncharged fluids. This first integral is subject 
to two integrability conditions expressed as nonlinear integral equations. We 
transform the integrability conditions, in Section 4, into a new system of dif- 
ferential equations which can be integrated in terms of quadratures. In Section 
5 we comprehensively investigate the nature of the factors of the quartic aris- 
ing in the quadrature. Finally, in Section 6, we discuss the results obtained 
and comment on some of the physical aspects. The approach of Khalique et al 
[T5] using a group classification may be helpful in providing further insights 
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in future work. 



2 Field equations 



We consider the shear-free motion of a spherically symmetric perfect fluid 
in the presence of the electromagnetic field. We choose a coordinate system 
x l = (t, r, 9, 0) which is both comoving and isotropic. In this coordinate system 
the metric can be written as 

dg 2 = _ e 2,(t,r) dt 2 + e 2X(t,r) ^2 + r 2 ^2 + sin 2 ^2)] ^ 

where v and A are the gravitational potentials. We are investigating the general 
case of a self-gravitating fluid in the presence of the electromagnetic field 
without placing arbitrary restrictions on the potentials. For this model the 
Einstein equations are supplemented with Maxwell equations. The Einstein 
field equations for a charged perfect fluid can be written as the system 



p=*£-^fa + >i + ^)-£* (2a) 

P = jl, (-3A? " 2A« + 2,A) + j- x (A, 2 + 2u r X r + ^ + ^) 

+^4T (2b) 



E 2 



r 4 e 4A 



(2c) 



= Ur A t - X tr (2d) 
Maxwell's equations yield 



E = r 2 e x - V $r, E r = ar 2 e 3X (3) 

In the above p is the energy density and p is the isotropic pressure which are 
measured relative to the four-velocity u a = (e~ u , 0, 0, 0). Subscripts refer to 
partial derivatives with respect to that variable. The quantity E = E(r) is 
an arbitrary constant of integration and a is the proper charge density of the 
fluid. We interpret E as the total charge contained within the sphere of radius 
r centred around the origin of the coordinate system. Note that $ r = F 10 is the 
only nonzero component of the electromagnetic field tensor F ab = <p b . a — <p a . b 
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where a = ($(£, r), 0, 0, 0). The Einstein-Maxwell system ([I])-© is a coupled 
system of equations in the variables p, p, E, a, v and A. 



The system of partial differential equations (j2J) can be simplified to produce 
an underlying nonlinear second order equation. Equation (12dl) can be written 

as 

v r = (lnA t ) r 

Then ([2b]) and (T2c]) imply 



e A [ A T 



+ 



2E 2 e 



2„-\ 







J t 



and the potential v has been eliminated. The Einstein field equations (j2J) can 
therefore be written in the equivalent form 



p = 3e 2h -e 



1 



e I X rr — A — 



A*e 3A 
e u = Xte~ h 



2A ( 2A rr + A 2 H 



E 2 



X xl + 



2X r 



— e 



3A+2h. 



2E 2 

r 4gA 



r 4gA 



(4a) 

(4b) 
(4c) 
(4d) 



for a charged relativistic fluid. In the above h = h(t) and F = F(r) are 
arbitrary constants of integration. Equation f!4dj) is the condition of pressure 
isotropy generalised to include the electric field. To find an exact solution of 
the field equations, we need to specify the functions h, F and E and solve 
equation (J4d|) for A. We can then compute the quantities p and p from (Jlal) 
and (jib"]) , and a follows from ©. 



It is possible to write ( I4dl) in a simpler form by eliminating the exponential fac- 
tor e A . We use the transformation, first introduced by Faulkes [17] for neutral 
fluids, which has the adapted form 



x 



y 



9(*) = ^ 



Then (]4dl) becomes 

y xx = f(x)y 2 + g(x)y 3 (5) 

which is the fundamental equation governing the behaviour of a shear-free 
charged fluid. Observe that (jSJ) is a nonlinear partial differential equation 
since y = y(t, x). When g = then y xx = f(x)y 2 for a neutral fluid which has 
been studied by Maharaj et al [19] and others. 
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3 A charged integral 



It would appear that we need to specify the functions f(x) and g(x) to inte- 
grate However it is possible ab inito to generate a first integral without 
choosing f(x) and g(x) if we generalise a technique first suggested by Srivas- 
tava [TT] . and extended by Maharaj et al [19]. The first integral generated 
is subject to a system of integral equations in f(x) and g(x) which can be 
rewritten as differential equations. 

Rather than choose f(x) and g(x) we seek general conditions that reduce the 
order of ([5]) to produce a first order differential equation. We can formally 
integrate ([5]) to obtain 



Vx = J f(x)y 2 dx + J g{x)y 3 dx 

= fiV 2 + 9iV 3 - 2 J fiyy x dx - 3 J giy 2 y x dx (6) 

For convenience we have used the notation 



f(x)dx = fr, J g(x)dx = g r 
Integrating fiyy x by parts and utilising ([5]) gives the result 

/ fiVVxdx = fnVVx - J fiivldx - J ff n y 3 dx - J gf n y A dx (7) 
We substitute (J7]) in to obtain 

Vx = fiy 2 + giy 3 - VnWx + 2 J fny 2 x dx 

+2 J ff n y 3 dx + 2 J gf n y 4 dx -3 J gjy 2 y x dx. (8) 



We continue this process and evaluate the integrals of fny 2 ., ffny and gfuy 
in ([8]) using integration by parts. Eventually we arrive at the expression 



yx = fry 2 + giy 3 - 2f n yy x + 2f in y 2 x + 2 y 3 + 2 {gfn), y 4 



2 -f I [2//jjj + 3 (//„), + l gi 



dx 



I /■ dx 



[gfrn + 2 (gfn);) \ dx (9) 
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For a meaningful result the integrals on the right hand side of (Q must be 
eliminated. 

We note that these integrals can be determined if 2f f ni + 3 (////) 7 + \gi and 
gfni + 2 (gfn)i are constants. This observation yields the following result 



r (t) = -y x + fry 2 + giy 3 - 2f n yy x + 2f IH y 2 x + 2[{ff II ) I - -K ]y 3 

+ [2(^ // ) / -ir 1 ]y 4 (10) 

subject to the integrability conditions 



K = 2ffr n + 3(ff II ) I + lgi (11a) 
K 1 =gf III + 2(gf II ) I (lib) 

where K and K\ are constants, and the quantity r (t) is an arbitrary function 
of integration. We have therefore established that a first integral of the field 
equation dSJ) is given by (fTO]) subject to conditions (|TT|) which are integral 
equations. On the surface it appears that the functions f(x) and g(x) are 
free. However equations (11 lap and (lllbj) effectively determine the forms of the 
functions f(x) and g(x); they are constrained by the integrability conditions 



4 Integrability conditions 



It is not easy to solve the nonlinear integral equations (fTTl) . However we can 
transform these equations into an equivalent system comprising a first order 
and a fourth order ordinary differential equation which are more convenient 
to work with. 

We let 

Zm = T 

so that /// = F x , fi = Fxx and / = J- XX x- Then it is possible to rewrite flllbl) 

as 

(gT) x + 2gT x = (12) 

Note that the integral equation (lllbl) has been transformed to a first order 
differential equation in T . Equation f|T2l) is integrable and we obtain 

g = /CoJ-- 3 (13) 

where T = Fi^x) and /C is an arbitrary constant. 
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Similarly we can eliminate g in (11 lap , with the help of ({TBI) , to get the result 

5 3 
TT XXXX + -T X T XXX = -^/Co^- 3 (14) 

Therefore the integral equation filial) has been transformed to a fourth order 
differential equation in T . Equation ( fT4l) can be integrated repeatedly to yield 



T- x =Kl + 1C 3 J F~V 2 dx + JC 2 (J F~ 3/2 dx^ 2 - i/d (J F'^dxX 

+ ^/C (/^ 3/2 rfx) 4 (15) 

where the /Q are arbitrary constants. 

We can rewrite (ITS"]) in a simpler form if we let 

J J=- Zj2 dx (16) 



so that 

Then we can write (HHP as 



u 



u. 



3/2 



/ 1 1 \ 3 / 2 

u x = [JC 4 + K 3 u + K 2 u 2 - -JC lU 3 + ^/C M 4 j (17) 

which is a first order equation in u. The equivalent integral representation is 

du 

x — x Q 



(£ 4 + /C 3 m + JC 2 u 2 - (l/6)/C lM 3 + (1/32)/C m 4 ) 3/2 

where x is a constant. The quadrature <HM can be evaluated in terms of 
elliptic integrals. We can summarise our result as follows: the first integral 
(TIP)), with g = ICqJ 7 ^ 3 , f = T xxx and T given by / TTgj) via /[T5\) . represents a 
particular class of solutions of (TJJ). 



To obtain solutions in closed form, satisfying the integrability conditions ( ITT]) , 
we need to evaluate the integral ( fl8l) . Particular solutions in terms of elemen- 
tary functions are admitted. In general the solution will be given in terms of 
special functions. We can express the solutions to ( TTTj) in the parametric form 
as follows 



f[x)=7 xxx (19a) 

g(x) = K Q F- 3 (19b) 

u x = ^ 2 = [G\u)}- 1 (19c) 

x-x = G(u) (19d) 
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where we have set 



G{u) = / — m (20) 

' (/C 4 + /C 3 m + lC 2 u 2 - (l/6)/dM 3 + (l/32)/C w 4 ) 3/ 



If we set g = (which forces /Co = 0) then the charge vanishes and the system 
( TT9"j) becomes 

f{x)=J c xxx (21a) 
Mx = jr-3/2 = ( 21b) 
x _a; = G( u ) (21c) 

where 



G(m) = / ^, (22) 

■/ (/C 4 + /C 3 m + /C 2 m 2 - (l/6)/C lM 3 ) 3/2 

This corresponds to the results found by Maharaj et al [19] for a neutral 

shear-free gravitating fluid. Thus their first integral is contained in our class 

of charged models (fT9|)-(l20|). 



5 Particular solutions 



Nine cases arise from the solution (TT9]) — (T2Q|) depending on the nature of the 
factors of the polynomial JC4 + JC311 + K2U 2 — (l/6)/Ciw 3 + (l/32)/C ^ 4 - 



5.1 Case I: One order-four linear factor 

If JC4 + IC3U + K2U 2 — (1/6)/Cim 3 + (l/32)/Co« 4 has one repeated linear factor 
then we have 

£ 4 + JC 3 u + K, 2 u 2 - (l/6)/C!M 3 + (1/32)/C m 4 = (a + bu)\ b^O 
We evaluate the integral in (1201) to obtain 



G(u) = -^-(a + bu)- 5 (23a) 
56 
24 

/(x) = -(56) 4 / 5 (x-x )- 11/5 (23b) 

(?(x)=/Co(56)- 12 / 5 (x-a;o)- 12/5 (23c) 

In this case it is possible to invert the integral (1TB]) and then write u = u(x). 
The first integral ( fTUl) has the form 
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r (t) = -y x - ^(56) 4/5 (x - x )- b/0 y 2 - ^/C (56) 



-6/5 2 



-12/5 



15 



7 



x — x J 



- 7/ v 



^(56) 4 / 5 (x - x r 1/5 yy x + y(56) 4 / 5 (x - x ) 4/5 ^ 
^(56) 8 / 5 (x - x r 7/5 - g/C Q (56) (x - x r 7/5 



-^/C (5fe)- 8 / 5 (x-x ) 



-8/5 4 

y 



(24) 



where we have used the functional forms in ( 123]) . The first integral ( 1241) corre- 
sponds to a shear-free spherically symmetric charged fluid which does not have 
an uncharged limit since /Co 7^ 0. If /Co = then the polynomial becomes cubic 
which is a contradiction. The charged integral (124"]) (£? 7^ 0, /C 7^ 0, 6 7^ 0) is 
a new solution to the Einstein-Maxwell field equations. 



5.2 Case II: One order-three linear factor 



If /C4 + K.3U + K2U 2 — (l/6)/Ciu 3 + (1/32)/C m 4 has two linear factors, one of 
which is not repeated, then we have 

£ 4 + IC 3 u + K, 2 u 2 - (l/6)/Ci<u 3 + (l/32)/C u 4 = (a + feu) (u + cf 



We use the computer package Mathematica [20] to determine the integral in 
(EUD to obtain 



G{u) 



2^J (a + bu)(u + c) 
35(a - 6c) 5 
136(a - 6c) 2 5(a - 6c)' 



356 4 9 36 3 2 96 2 (a - 6c) 



a + bu u + c 



{u + cf 



(u + cY 



{u + c) 4 



(25) 



expressed completely in terms of elementary functions. In this case, if g = 0, 
/C = and 6 = 0, then (1251) becomes 



-7/2 



(26) 



GW = «- 3/2 (ti + c) 

and hence using (l2Tj) we find 

^)=a 2/7 (^)(4) 6/7 (x-xo)- 5/7 (27) 
Note that (1271) is related to the result obtained by Maharaj et al [T9] . 
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Again setting g = 0, K\ = 0, in (fit)]) we get 



Mt) = -y* + fiv 2 + aiy 3 - Vim, + Vinyl + Wfnh - 



which was the first integral for uncharged matter found by Maharaj et al [19] . 
Also observe that if g = 0, ifi = 0, f(x) = (ax + b)~ 15/7 then (03JD yields 



which was found by Srivastava [TT] . Also with # = 0, ifi = 0, /(x) = x 15 ^ 7 in 
(fTOj) (or if we set a = 1, 6 = in (J2HD) we have 



which was established by Stephani [2T]. Therefore the first integral (1101) is a 
charged generalisation of the particular Maharaj et al [19J, Srivastava [llj 
and Stephani [21] neutral models. 



5.3 Case III: One order-two linear factor; one order-one quadratic factor 

If /C4 + /C3U + K.2U 2 — (1/6)/Cim 3 + (l/32)/Co« 4 has two factors, one linear and 
repeated and the other is irreducible to linear factors, then we have 




(28) 




/C 4 + IC 3 u + K 2 u 2 - (1/6)/Cim 3 + (l/32)/C « 4 




The function (l2"0l is integrated to obtain 
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G{u) 



f 1 5(6 - 2cd) 

\(a — bd + cd?) v? 2(a — bd + cd?)u 

15(6 - 2cdf - 62c(6 - 2cd) 2 (a - bd + cd 2 ) + 2Ac 2 {a - bd + cd 2 ) 2 

2(a - 6d + cd 2 ) [4c(a - bd + cd 2 ) - (6 - 2cd) 2 } 
c(6 - 2cd) [15(6 - 2ccf) 2 - 52c(a - bd + cd 2 )] u ' 



2(a-bd + cd 2 )A 
1 



x 



2yJ(a -bd + cd 2 ) + (6 - 2cd)u + cu 2 

15(6 - 2cd) 2 - \2c{a -bd + cd 2 ) 
+ 8(a - bd + cd?f X 

du 



f , (29) 
J uJ a — bd + cd 2 + (6 — 2cd)u + cu 2 



where A = 4(a — bd + cd 2 )c — (6 — 2c<i) 2 and the integral on the right hand 
side can be expressed in terms of elementary functions. The exact form of 
the integral depends on the signs of a — bd + cd 2 and A (see Gradshteyn and 
Ryzhik [22], equations 2.266 and 2.269.6). 



5.4 Case IV: One order-two linear factor; two order-one linear factors 

With one repeated and two non-repeated linear factors we have 

£4 + JC 3 u + }C 2 u 2 - (1/6)/Cim 3 + (1/32)/C m 4 = (a + 6m) (cu + d)(u + e) 2 

In this case the expression for the integral in (|20~!) can be evaluated with 
the help of the computer package Mathematica [20] . The resulting expression 
is expressible in terms of only elementary functions. This expression is very 
lengthy and not illuminating, and is therefore not included in this work. 

5.5 Case V: Two order-two linear factors 

If /C4 + /C3M + K2U 2 — (1/6)/Cim 3 + ( 1/32) /Co w 4 has two linear factors each of 
which is repeated, then we have 

/C 4 + /C 3 m + }C 2 u 2 - (l/6)/Ci« 3 + (1/32)/C m 4 = (a + bu) 2 (u + c) 2 . 

The integral in (1201) may be easily determined so that 
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G(u) 



Qb 2 In 



u + c 3b 2 (a-bc) b 2 (a-bc) 2 3b(a - be) 



a + bu 



+ 



a + bu 



+ 



2(a + 



+ 



u + c 



(a - 6c) 2 
2(w + c) 2 



(a - 6c) 5 



(30) 



Thus for the case of two order-two linear factors the integral can be expressed 
completely in terms of elementary functions. 



5.6 Case VI: No repeated linear factors 



If /C 4 + IC 3 u + K 2 u 2 — (l/6)/CiM 3 + (1/32)/C m 4 has no repeated linear factors, 
then we have 

/C 4 + /C 3 m + /C 2 m 2 -(1/6)/CiM 3 + (1/32)/CoM 4 = e(a+u)(b+u)(c+u)(d+u),e^ 

In this case we obtain, in terms of elementary functions and elliptic integrals, 
the result 



G{u) 



2e" 3 / 2 



(a - b)^J(a + u)(b + u)(c + u)(d + u) L ( 6 ~ c )(° ~ d ) 



(a + u)(b + u) 



(b - dy 



+ 



(b — d)(c — d) (a — c)(c — d) + 
1 1 



(b-c)(b-d) (a-d)(b-d) 
4e -3/2 



b + u 
a — c 
1 



2(d + u) 



(a — b)(a — df 



{a-b)Vb^d 



(b-c)(b-d) 



+ 



y/a — c 



(a — d) 2 (c — d)^a — c (a — b)(b — c) 2 (b — d) 



+ 



a — b — c + d 



E(a,p) 



{c-df{a~cf/ 2 {b-c) 
2e- 3 / 2 

+ {a- c y/ 2 {b-dy/ 2 {b-c){a-d) x 

2(a + b-c-d) 2 (a-b-c + d) 2 
{b — c) (a — d) (a — b) (c — d) 

where we have let 



F(a,p), (0<d<c<b<a) (31) 



a = arcsm 



(a — c)(d + u) 
\ (a — d){c + u) 



P 



{b — c)(a — d) 
(a — c)(b — d) 



In (!3T|) . F(a, p) is the elliptic integral of the first kind and E(a, p) is the elliptic 
integral of the second kind. This result is similar to one of the results obtained 
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by Maharaj et al [19]. However their uncharged model is not regainable from 
the expression above as the polynomial here is necessarily quartic. 



5. 7 Case VII: One order-two quadratic factor 



If /C 4 + IC 3 u + IC 2 u 2 — (l/6)/CiM 3 + (1/32)/C m 4 has one repeated quadratic 
irreducible factor, then we have 

/C 4 + /C 3 m + K, 2 u 2 - (l/6)/C!U 3 + (1/32)/C m 4 = (a + bu + cu 2 ) 2 

In this case we obtain 



G{u) 



b + 2cu 



Aac 
+ 



-b 2 

6c 2 



2(a 



1 

bu + cu 2 ) 2 
du 



+ 



3c 

(4ac — 6 2 )(a + bu + cu 2 



(4ac — 6 2 ) 2 7 a + 6m + cu 2 



which can be expressed in terms of only elementary functions. The exact form 
of the integral depends on the sign of 4ac — b 2 (see Gradshteyn and Ryzhik 
[22], equations 2.172 and 2.173.2). 



5.8 Case VIII: Two order-one quadratic factors 



With two non-repeated quadratic factors we have 

K A + /C 3 m + K 2 u 2 - (l/6)/Ci« 3 + (1/32)/C m 4 = (a + bu + cu 2 )(d + eu + u 2 ) 

In this case the expression for the integral in ( 1201 , using the computer package 
Mathematica [20], is obtainable but is not included in this work as it is very 
lengthy. It may be expressed in terms of elementary functions and elliptic 
integrals. 



5.9 Case IX: One order-one cubic factor 



If /C 4 + /C 3 m + /C 2 m 2 — (l/6)/CiM 3 + (1/32)/C m 4 has one irreducible cubic factor, 
then we have 

/C 4 + JC 3 u + IC 2 u 2 - (l/6)/Ciu 3 + (1/32)/C m 4 = (a + bu + cu 2 + du 3 )(e + u) 

The integral in (1201) can again be found with the help of the computer package 
Mathematica [2U] . It is given in terms of elementary functions, elliptic integrals 
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and special functions. However it is so lengthy that it is also not included in 
this work. 



6 Discussion 

In this paper we have modelled the behaviour of shear-free charged fluids, 
and reduced the solution of the Einstein-Maxwell system of field equations to 
a single nonlinear partial differential equation. By treating this equation as 
an ordinary differential equation, a first integral was found using elementary 
methods. It is remarkable to note that the first integral is obtainable without 
specifying the arbitrary functions contained in the governing equation. The 
first integral is subject to a system of two integral equations which were re- 
placed by a system of two differential equations which can be integrated up to 
a quadrature. Consequently we have found a new class of parametric solutions 
to the Einstein-Maxwell system for a charged gravitating shear-free fluid. The 
new solution is given by the parametric equations (}T9"]) - (j2"u]) . 

A detailed analysis of the factors of the quartic arising in the quadrature 
was performed. Two cases of interest arise. Firstly, we are in a position to 
explicitly invert the quadrature when there is one repeated linear factor and 
explicitly write the first integral. Then the model has to be necessarily charged. 
We believe that this is a new result. Secondly, we can explicitly invert the 
quadrature when there is one order-three linear factor. This case contains 
that of vanishing charge and we regain the results of Maharaj et al |19j . 
Srivastava [11] and Stephani [21]. In the remaining cases the functions G(u) is 
a complicated combination of elementary functions or/and special functions. 
In these remaining cases it is not possible, except maybe for special parameter 
values, to invert the integral and write expressions for f(x) and g(x) explicitly. 

We make certain points, related to the mathematics, to clarify the approach 
followed in this paper. Firstly, the Einstein-Maxwell system has been studied 
extensively in the past with the objective of finding exact solutions. We have 
considered earlier treatments and, in particular the comprehensive analyses 
of Krasinski [13J and Stephani [2T], and have not found any reference to the 
first integrals established in this paper. Secondly, we have generated the first 
integral ffTUl) mathematically following the approach of Maharaj et al [19J for 
uncharged fluids. We believe that this is an elegant approach and may be used 
in other investigations for the gravitational field or other physical systems. 
Our new class of solutions may be useful in this context, and could provide a 
deeper insight into the behaviour of the gravitational field. A comprehensive 
mathematical analysis of the integrability properties of ([5]) using the symmetry 
properties of the equation may provide further solutions and insights. For 
example the treatment of Halburd [23], for the uncharged shear- free case, 
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established an equivalence with the generalised Chazy equation and provided 
a new class of integrable equations. 



We now make certain comments, related to the physics, related to our results. 
Firstly, a natural question is whether the models generated here are physi- 
cally meaningful. It is difficult to perform a general qualitative study of the 
physical features of the models because of the complexity of the functions 
involved. However we have found explicit charged first integrals and earlier 
uncharged first integrals are regained. Other parameter values may also lead 
to acceptable models. These factors point to physical reasonability. Secondly, 
the line element ([I]) is written in terms of isotropic and comoving coordinates; 
the spacetime is shear-free. The introduction of charge, and assumptions made 
in the integration, do not affect ([I]) so that the spacetime remains shear-free. 
Thirdly, charged shear-free models in the presence of heat flow are of crucial 
importance in relativistic astrophysics |7;8j and influence the range of temper- 
ature profiles of the models discussed here. In [25] it was shown that, for the 
most general spherically symmetric line element with acceleration, the causal 
transport equation reduces to 

0{qB)jr- + A(qB) = _ a ^^k (32 ) 



where 



K = 7 T 3 r c r = ( ^ ] r c = (3T~ a (33) 



a 



with 7>0,a>0,/3>0 and a > 0. Then it is possible to integrate this 
equation for general metric functions A and B. When (3 = 0, all noncausal 
solutions of (|32|) are given by 



M T )4-a = Z — 1 f A 4 -*qB 2 dr + F(t) a^A (34) 
a J 

\n(AT) = — [ qB 2 dr + Fit) a = 4 (35) 
a J 

where F(t) is an arbitrary function of integration which is fixed by the ex- 
pression for the temperature of the star at its surface S. For causal solutions 
((3 7^ 0) two solutions are provided. In the case of constant mean collision 
time, ie. a = 0, ( !32l is simply integrated to yield 



(AT) 4 = -- 
a 



f3 J A 3 B(qB), t dr + J A 4 qB 2 dr 



+ F{t) (36) 



The only nonconstant mean collision time solution is given for a = 4: 
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(AT) 4 = exp (- J ^dr^j J A 3 B(qB), t exp ^ J ^dr^j dr 

+ F(t)exp(-J ^dr^j (37) 

Our results can be incorporated into the above framework with A = e u and 
B = e x in the astrophysical context. 

The main objective of this paper was to show that an earlier method to obtain 
uncharged first integrals is extendible to the Einstein-Maxwell system of equa- 
tions. We have shown that this is possible and particular charged first integrals 
have been explicitly obtained. It would be beneficial to study the models gen- 
erated in terms of the original metric, at least in particular cases for simple 
forms of f(x) and g(x), and to evaluate the evolution of the charged relativis- 
ts fluid in terms of the fluid and electromagnetic variables. Other quantities 
requiring attention are the spacetime structure, energy conditions, Riemann 
invariants and causality. This is outside the scope of our present treatment 
and will be considered in future work. 
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